On the interpretation of echelle diagrams for solar-like oscillations 

Effect of centrifugal distortion 



J.C. Suarez^'^ and M.J. GoupiP and D.R. Reese^ and R. Samadi^ and F. Lignieres^'^ and 

M. Rieutord'^'^ and J. Lochard^ 



ABSTRACT 

This work aims at determining the impact of slow to moderate rotation on the regular patterns 
often present in solar-like oscillation spectra, i.e. the frequency spacings. We focus on the 
well-known asteroseismic diagnostic echelle diagrams, examining how rotation may modify the 
estimates of the large and small spacings, as well as the identification of modes. We illustrate 
the work with a real case: the solar-like star 77 Bootis. 

We study a main sequence 1.3Mq star as a typical case. The modeling takes into account 
rotation effects on the equilibrium models through an effective gravity and on the oscillation 
frequencies through both perturbative and non-perturbative calculations. We compare the results 
of both type of calculations in the context of the regular spacings (like the small spacings and the 
scaled small spacings) and echelle diagrams. We show that for echelle diagrams the perturbative 
approach remains valid for rotational velocities up to 40 — 50kms^^. 

We show that for the rotational velocities measured in solar-like stars, i.e. vsini up to 20-30 
kms~^, rotation effects must be taken into account in the modeling for a correct interpretation of 
the observed oscillations. In particular, theoretical oscillation frequencies must be corrected up 
to the second-order in terms of rotation rate, including near degeneracy effects. For rotational 
velocities of about 16 kms~^ and higher, diagnostics on large spacings and on modal identification 
through echelle diagrams can be significantly altered by the presence of the m components 
of the rotationally split modes. We found these effects to be detectable in the observed frequency 
range. Analysis of the effects of rotation on small spacings and scaled small spacings reveals that 
these can be of the order of, or even larger than surface effects, typically turbulence, microscopic 
diffusion, etc. Furthermore, we show that scaled spacings are significantly affected by stellar 
distortion even for small stellar rotational velocities (from 10-15 kms^^) and therefore some care 
must be taken when using them as indicators for probing deep stellar interiors. 

Subject headtngs: stars: evolution -stars: individual(7; Bootis) — stars: interiors - stars: oscillations 
(including pulsations) stars: rotation -stars: solar-type 

1. Introduction 

Institute de Astrofi'sica dc Andalucia (CSIC), Granada, 

Spain. Solar-like oscillations are detected in stars with 

^Observatoire de Paris, LESIA, CNRS UMR 8109, Outer convective zones, that is. in s tars similar 

92195 Meudon, France. to the Sun (see iMichel et ahl boost and ref- 



3Universite de Toulouse, UPS, Laboratoire erences therein) as w ell as in red giant stars 



d'Astrophysique de Toulouse- Tarbe (LATT), 31400 VTZ — t-,. , . ~ — Tl U„„rJ\ j ■ • n 

Toulouse France. (IDe Ridder et al.l I2009I I and eve n m massive B- 



4CNRS, Laboratoire d'Astrophysique de Toulouse- tyM stars as recently shown by [ Belkacem et al. | 

Tarbes (LATT), 31400 Toulouse, France. (|2009l ). These osciUations are due to intrinsically 



stable pressure (p) modes which are excited by 
turbulent motions in the convective outer layers. 
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This stochastic process excites several modes in a 
star over a broad frequency range weh above the 
fundamental radial mode frequency in a regime 
where frequencies show regular patterns in the 
power spectrum. This gives access to an infor- 
mation content based on different seismic signa- 
tures which may constrain fundamental stellar 
properties (namely mass, age, and radius) and 
probe the internal stellar structure to unprece- 
dented accuracy (see, for instance. lGough| | l987 



Cunha fc Metcalfell2007l : iTeixeira et al.ll2009i r 



On the other hand, the uninterrupted long ob- 
servational runs required together with the ex- 
pected small oscillation amplitudes has made the 
detection of solar-type stars technically challeng- 
ing targets in particular, for asteroseismic space 
missions. A wealth of data is already a vailable 
from MQgT (iMatthewsl Il998l). CoRoT (iBaglin 



20031: iMichel et al.ll2008t lAppourchaux et al.ll2008t 



Belkacem etal]|2009l ) and Kepler (?). 

Low mass stars are generally slow rotators (in 
most cases {vsini) < 20kms^^) and the influ- 
ence of rotation on the oscillation frequencies is 
globally small. For slow rotators, the geometri- 
cal structure is represented by a single spherical 
harmonics Yp{9, (/>), and the frequencies by h'n,£,m 
with n being the radial order of the mode. Ro- 
tation lifts the degeneracy causing the modes to 
split into m 7^ multiplets. The m ^ are 
usually computed using a first-order perturbation 
formalism. Rotation also affects, indirectly, the 
m = mode frequencies as it modifies the inter- 
nal structure through rotationally induced mixing 
(iGoupil fc Talonll2002l: iMathis et al.ll2007l ) . In the 
same line, Carrier et al.l ( 2005 ) studied the solar- 
like star r] Bootis assuming a surface rotational 
velocity of about Wrot = 3 — 7kms^^. The authors 
found that the effects of including rotationally- 
induced mixing was small and the rotation is not 
strong enough to impose a rigid rotation. In par- 
ticular, they found that the core in this star rotates 
three times faster than the outer layers. 

Distortion due to the centrifugal force can have 
a much larger impact on the o scillation freq uen- 
cies even for slow rotators (see Goupil 20091 for 
a review on rotation effects on the oscillation fre- 
quencies). Such an effect is stronger for p modes 
with small inertia. These modes propagate mainly 
through the outer layers of the star. Therefore, 
their frequencies are more sensitive to changes in 



the surface physical properties, where the centrifu- 
gal force becomes more efficient. The present work 
examines and evaluates the direct effects of rota- 
tion on the frequencies in this high-frequency do- 
main, and its impact on the different seismic diag- 
nostic techniques based on the asymptotic proper- 
ties of the oscillations. More specifically, we ana- 
lyze the case of different modes with close frequen- 
cies, such as £ — and £ = 2, which are system- 
atically near degenerate. Since the small spacings 
are affected by near degeneracy, echelle diagrams, 
often used as a seismic diagnostic technique for 
solar-like stars, should also be affected. 

We then first analyzed the effects of rotation on 
regular patterns observed in solar-like spectra, fo- 
cusing on the small and large spacings, which are 
the most relevant quantities for echelle spectro- 
grams. Even if solar-like stars are often slow rota- 
tors, their oscillation spectra show high-order fre- 
quencies, where the validity of perturbation theory 
for computing oscillation frequencies in absence of 
rotation might fail. This point was thus checked 
by comparing the perturbative calculations with 
recent non-perturbative calculations where the ef- 
fect of rotati on on the oscillat ions is fully taken 
into account ( Reese et al. 20061) . Taking all these 
previous steps into account, we investigated the 
effect of rotation on echelle diagrams. 

The paper is organized as follows: in Section [2] 
for sake of notation, we briefly recall the frequency 
spacings and seismic diagnostics generally used for 
asteroseismic analysis of solar-like stars which will 
be considered in the following sections. Then, the 
equilibrium stellar models and their correspond- 
ing oscillation computations are described in Sec- 
tion [21 The validity of a perturbative modeling 
of the effects of rotation on oscillations in solar- 
like stars is discussed in Section S) Section [5] de- 
scribes the effects of rotation for seismic diagnos- 
tics based on echelle diagrams. A discussion of 
the results, including an illustration based on the 
solar-like star -q Bootis and conclusions are given 
in Section [6] 

2. Seismic diagnostics for internal struc- 
ture 

In the case of solar-like oscillations, modes with 
high radial order n are detected. In this regime, 
regular patterns are observed in the power spec- 
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trum which are related to internal properties of 
the star. Analysis of such characteristics allows us 
to infer information about various physical pro- 
cesses, particularly those that modify the sound 
speed and/or the Brunt- Vaisala frequency. Sev- 
eral frequency combinations are then used to an- 
alyze high-order p modes: large and small fre- 
quency spacings and their corresponding scaled 
counterparts, second-order frequency differences, 
and even high-order freq uency differences (for a 
review on this topic see iGoupil fc Dupreti 12007 , 
and references therein). 

2.1. Frequency spacings 

The large (frequency) separation is hereafter 
denoted by Aj/„,f and defined by 



Aj/„ 



(1) 



The large spacing is sensitive to surface layers 
where Cs is small, and hence oscillation modes 
spend more time there compared to the inner re- 
gions of the star. This spacing is generally studied 
in order to detect sharp variations of the sound 
speed, such as those at the base of the solar con- 
vection zone, or those resulting from the presence 
of ionization regions. The mean large spacing 
over several modes in the asymptotic regime scales 
with the dynamical time scale (i?^/G'M)^/^ and is 
one of the pieces of information sought when con- 
structing echelle diagrams. 

The small (frequency) spacing, hereafter de- 
noted by (5i^„_£, is defined as follows 



(2) 



This spacing is sensitive to evolutionary effects 
(the age of the star for a given metallicity). Av- 
erage large (Af„_£) and small {5vnj) spacings 
(over radial order n) are used to estimate masses 
and ages of solar-like s tars a s first suggested by 
Christensen-Dalseaar dl (|l998h . 

In addition to the aforementioned frequency 
spacings, other frequency combinations are also 
used as efficient asteroseismic diagnostics: 



^0,l("-) = I'nfi - -j^{l^n-lA+Vn,l) (3) 
5lfl{n) = ]^{Vnfl + y-n+lfi) - l^n.l, (4) 



together with their corresponding scaled expres- 
sions 



n-l-1,0 



?'i,o(n-) 



'^i,o(?t) 



(5) 



which are sensitive to inner stellar zones (iRoxburgh fc Vorontsov 



2001 



Mazumdar fc Antiall2001 : Roxburgh h Vorontsov 



20031 see). 



2.2. The Echelle Diagrams for asteroseis- 
mic diagnostics 

One of the most frequently used techniques for 
seismic diagnostics of solar-like stars is the rep- 
resentation of the oscillation frequencies in the 
so-called echelle diagram (hereafter, ED). With 
such diagrams, one is able to identify the degree 
of modes and extract some information from the 
asymptotic properties mentioned above. 

Echelle diagrams are constructed by cutting the 
oscillation spectrum into frequency slices of size 
(Ai/n ^), and stacking them up. Then the ED con- 
sists in depicting the oscillation frequencies as a 
function of the same frequencies modulo {Avn,ii- 
In a first-order asymptotic regime, EDs should be 
composed of one vertical ridge per ^, if modes are 
m-degenerate. When rotation is considered, each 
vertical ridge corresponds to a given (^, m) mode. 

3. The seismic modeling 

The present work has been undertaken with the 
use of two types of models: 1D-2D polytropes and 
ID stellar models, the purpose and description of 
which are described in the following. 

3.1. Polytropes 

Polytropic models are used in this work for 
the comparison between perturbative and non- 
perturbative approaches for the oscillation com- 
putations (section H]). This is so because, to date, 
no realistic (properly deformed by rotation) mod- 
els of solar-like stars are yet available. 

A polytropic model of index n — i with AI = 
1.3 Mq and R — 1.276 R© was considered as the 
reference model. From that model, ID and 2D 
polytropic models with rotational velocities rang- 
ing from 8kms~ ^ to 48kms^^ were built (see 



iReese et al. 2006 . for more details on 2D poly- 
tropes). The characteristics of all these models 
are listed in Table [T] 
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Table 1 

Characteristics of the 1.30 models considered in this work (with the exception of Mr, 

WHOSE MASS IS 1.70 Mq). Pi, AND Vj ARE ID AND 2D POLYTROPIC MODELS, RESPECTIVELY. StELLAR 
models are denoted by Ui AND Si, FOR UNIFORM (CONSTANT ROTATION PROFILE AT GIVEN 
EVOLUTIONARY STAGe) AND SHELLULAR ROTATION (ROTATION PROFILE GIVEN BY EQ. [7|) MODELS, 
RESPECTIVELY (DETAILS IN SECTION 13. 3p . MODEL SUBSCRIPTS i MAKE REFERENCE TO THE ROUND 
VALUE OF THE SURFACE ROTATIONAL VELOCITY IN KMS~^ (fOR Mr, THE ROTATIONAL VELOCITY IS 

30kms^^ approximately). The radius of all the models is i? = 1.276 ± 0.005 R©. Columns 

REPRESENT RESPECTIVELY, FROM LEFT TO RIGHT, THE AVERAGE LARGE SPACING, THE ACOUSTIC 
FREQUENCY U^c (EO. [TT]) . THE VALUE OF t'max(EQ. [TO]) , THE ROTATION FREQUENCIES IN BOTH THE 
SURFACE AND THE CENTER OF THE MODEL, AND THE SMALL NUMBERS £ AND fj, (EQ. [9|) THE LATTER 
BEING CALCULATED AT THE VainK FREQUENCY OF EACH MODEL. 



Model 


/iHz 


mHz 


niHz 


fiRz 


fiRz 


e 

xlO"^ 


/i(l^max) 
XlO"^ 


^16 

" 32 
■pt 

48 


98.40 
98.40 
98.40 
98.40 
98.40 
98.40 







1.43 
2.92 
5.84 
5.84 
8.75 



1.43 
2.92 
5.84 
5.84 
8.75 




I. 85 
3.70 
7.40 
9.26 

II. 11 




Ps 

Pl6 
P32 
P40 

PiS 


98.40 
98.40 
98.40 
98.40 
98.40 






1.43 
2.92 
5.84 
7.24 
8.75 


1.43 
2.92 
5.84 
7.24 
8.75 


I. 81 
3.70 
7.40 
9.16 

II. 08 




Uqs 
Uw 
U32 
Um 

U4S 


108.60 
108.62 
108.65 
108.67 
108.70 


3.79 
3.79 
3.79 
3.79 
3.80 


2.25 
2.23 
2.24 
2.24 
2.24 


1.47 
2.99 
5.75 
7.20 
8.58 


1.47 
2.99 
5.75 
7.20 

8.58 


1.86 
3.79 
7.28 
9.12 
10.84 


0.66 
1.34 
2.58 
3.23 
3.83 


58 
5l6 

532 

540 

548 


108.61 
108.64 
108.66 
108.69 
108.72 


3.79 
3.79 
3.79 
3.79 
3.80 


2.24 
2.23 
2.24 
2.24 
2.24 


1.47 
2.92 
5.49 
7.29 
8.73 


I. 89 
3.79 
7.32 
9.69 

II. 48 


I. 86 
3.69 
6.91 
9.24 

II. 03 


0.66 
1.30 
2.45 
3.27 
3.90 


Mo 


109.12 


3.79 


2.25 


0.00 


0.00 








Mr, 


39.78 


1.07 


0.63 


2.54 


2.54 


9.11 


4.02 
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3.2. Stellar models 

Equilibriu m models were computed with th e 
CESAM code (|More]||l997 l :lMorel fc LebretoulbOOSl ). 
Our stellar models include the spherically av- 
eraged contribution of the centrifugal accel- 
eration through an effective gravitational ac- 
celerati on in the hydrostatic equati on accord- 
ing to Kippenhahn fc Weigert ( 199Cl() . Specifi- 
cally, the effective gravity is obtained as ges = 
g — Ac{r), where g is the local gravity, r the 
radius, and Ac{r) — |rf2^(r) the centrifugal 
acceleration (similar models have been previ- 
ously used, for instance, by [C oupil ct al. 200^ 



cies are close to each other. More specifically, in 
order to remain in the hypothesis of the pertur- 
bative approach, we assume the loose condition 
that near degeneracy effects must be corrected for 
modes with frequencies satisfying: 



\^nlm ^n'l'ra)\ 



(6) 



In addition, the frequency computation takes into 
account the presence of a radial differential rota- 
tion profile of the form 



r!(r) = n 1 + 770 (r) 



(7) 



Daszvnska-Daszkiewicz et al.l2002l : IChristensen-Dalsgaard fc ThoHipson 
1TT^nTT^~^~~nTrr7n 7777 \ n where si repre 



19991: ISuarez et all 120061 ) . Although the non 



spherical components of the centrifugal acceler- 
ation are not considered in the equilibrium mod- 
els, they are included in the oscillation compu- 
tations by means of a linear perturbation ana l 



ysis described iii iDziembowski fc Goodd ( 1992[ ). 



Soufi et all (|l998l ).and [Suarez et alj (|2006[) 



We consider two hypothesis about the inter- 
nal rotation profile: uniform rotation (solid body 
rotation), and shellular rotation (radial differen- 
tial rotation) with local conservation of the angu- 
lar momentum during the evolut ion of the stel- 
lar model (see ISuarez et al. 20061 for more de- 
tails) . The characteristics of both uniformly- and 
differentially-rotating models (W,; and Si, respec- 
tively) are listed in Table [T] We consider rota- 
tional velocities for the models spanning the range 
of V sin i observed in solar-like stars. 

3.3. Oscillations computation 

Oscillations were computed using the non- 
perturbative and perturbative approaches. The 
former consists in a numerical approach based 
on expansions of the equilibrium and oscillation 
variables on spherical harmonics for the angular 
dependence, and on Chebyschev polynomials for 
the radial dependence. Such computations were 
done for the 2D polytropic models described in 
Section 13.11 

For the perturbative approach, the adiabatic 
oscillation code FILOU was used ( Snared 2002 : 
ISuarez fc Goi^lioosl) . This code corrects the 
oscillations frequencies up to the second-order ef- 
fects of rotation. These include near degeneracy 
effects, which occur when two or more frequen- 



ere si represents the angular rotational veloc- 
ity at the surface and rjQ^r) a radial function rep- 
resenting shellular rotation. This profile is only 
considered in the radiative zone, in which the ro- 
tation rate decreases approximately with a power 
of r, whereas instantaneous transport of the angu- 
lar momentum is assumed in the convective core, 
implying thus a uniform rotation profile in that 
part of the star (an illustration of such a profile 
can be found in Suarez et al.ll2009[ ) . 



4. Validity of the perturbative approach 
for oscillation computations 

The perturbative approach is considered valid 
when the parameters 



= f]/(GAf/i?3)i/2 



(8) 
(9) 



are small, that is, when (1) the stellar structure is 
not significantly deformed by the centrifugal force, 
and (2) oscillation frequencies are much larger 
than the angular rotation rate, respectively. More 
specifically, the stellar deformation, which scales 
as e^/^, needs to be small in comparison to the 
mode wavelength, which scales as l/vn^g. This im- 
plies that the use of a perturbative approach may 
fail when increasing the rotational velocity of the 
star, and this failure is expected to come first for 
high radial-ord er frequency modes, i . e. high fre- 
quen cy modes ( Lignieres et al. 2006t Reese et al 



20061) . 



It is thus relevant to investigate the validity of 
the perturbative approach for the oscillation com- 
putations, as well as to identify the possible effects 
that it may provoke on the echelle diagrams. To 
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do so, it is necessary to compare EDs constructed 
from synthetic oscillations computed with both 
approaches. In order to try to remain within the 
nominal limits of validity of perturbation theory, 
models were built with small values of e and /i, 
around 10~^ and smaller (see Table [T]). 

For the sake of shortness and simplicity, P 
and NP will indicate the perturbative and non- 
perturbative cases, respectively. For instance, 
echelle diagrams computed using the perturba- 
tive approach will be denoted by PEDs (Pertur- 
bative Echelle Diagrams), and those computed 
with a non-perturbative approach will be denoted 
by NPEDs (Non- Perturbative Echelle Diagrams). 
For both P and NP models, uniform rotation is 
assumed. 

The comparison between both types of dia- 
grams including all the m components of modes 
shows that these are very similar for rotational 
velocities lower than 32kms~^ approximately (see 
Figure [Ij. For higher rotational velocities (e > 
0.09), the P results should be interpreted with cau- 
tion. 

Nevertheless only slight differences (for this 
scaling, i.e. a few /xHz) between PED and NPEDs 
come out for f = and £ — 1, to = modes in the 
high-frequency region {i'n,i > 2mHz). This means 
that, at the scale of echelle diagrams, both P and 
NP can be considered as quasi-identical. This is 
also the case when comparing the P and NP large 
frequency spacing. 

For the small spacings, we find very similar re- 
sults from both approaches, as can be deduced 
from the comparison between the top and bottom 
panels in Figure [5] For (5j^o,2, P spacings show 
the same behavior as a function of the frequency 
as their NP counterparts. Indeed, the results are 
quasi-identical for models with rotational veloci- 
ties up to 40kms^^. Deviations of P with respect 
to NP spacings are more evident for 48kms~^, 
showing differences of around 1 fiRz, which are vis- 
ible for frequencies higher than 1.5 mHz and which 
increase with the oscillation frequency. Similarly, 
the analysis of P and NP results for Si'i q are sim- 
ilar, with increasing differences from 1.5 mHz to 
higher frequencies. 

Furthermore, when near degeneracy effects are 
not considered, the resulting small spacings devi- 
ate from the NP results even for the lowest ro- 



tational velocities considered here. This is shown 
in Figure [2] (middle panels) , where it can be seen 
that such deviations increase with the frequency 
and the rotational velocity. This can be explained 
by the dependency of the near-degeneracy correct- 
ing terms upon, among other terms, the frequency, 
the rotation profile, and structure terms a ffected 
by the centrifugal force (jSuarez et al.ll2006l ). This 
issue, together with a more detailed and quantita- 
tive comparison, is currently in progress, and the 
results will be publish ed in a forthcoming paper 
( Ouazzani et al. 2010l in prep.). 

Several important consequences can be ex- 
tracted from these results: (1) It is correct to use 
the P approach for studying the effect of rotation 
on EDs (up to e = 0.09), (2) even for low rotational 
velocities, second-order effects of rotation, includ- 
ing near degeneracy effects, must be included in 
the oscillation computations, and (3) there are no 
artificial features in the PEDs which could jeop- 
ardize the correct physical interpretation of the 
seismic diagnostics performed with them (even at 
a small scale) as long as near degeneracy effects 
are properly taken into account. 

5. EDs of stellar models 

Following the conclusions on the validity of the 
perturbation methods given in the previous sec- 
tion, we computed the oscillations of the Ui and 
Si models (Table [T]) including near degeneracy ef- 
fects. It can be shown that EDs of Si and Ui mod- 
els are very similar, so for simplicity, discussion 
will be focused on Ui models. For brevity, ridges 
in the EDs are specified with the notation {£,m). 

Figure [3] shows the echelle diagrams calculated 
for the three Ui models of Tabled] At first sight, 
one notices that second-order effects of rotation 
do not change substantially the overall shape of 
the EDs, which remains dominated by the internal 
stratification of the star and the presence of hydro- 
gen and helium ionization zones. However, signif- 
icant variations (within the global ED structure) 
are found for increasing rotation rates 51, and for 
different frequency domains. In fact, such varia- 
tions should be considered in terms of (e, n) values, 
rather than a function of the rotational velocity. 
It can be shown, for instance, that for similar ro- 
tational velocities but different masses, EDs show 
different results. Here, for clarity, the discussion 
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is nevertheless held with respect to Wj-ot- 

For Uqs, the multiplets with different £ remain 
quite distant from each other. In particular the 
frequency gap between the £ = and 2 multiplets 
is around 20 /iHz, and between the £=1,3 multi- 
plets around 10 /iHz. This roughly corresponds to 
the spacing between to = ridges for those multi- 
plets. For such a low rotational velocity, first-order 
effects of rotation, which are proportional to mQ, 
dominate, thereby explaining the symmetry of the 
multiplets. In this case, the rotational splitting is 
1.47 /zHz (see Table [1]), which represents approx- 
imately the 7% and 14% of the average distance 
between the 0-2 and 1-3 multiplets, respectively, 
and even more importantly, about 1% of the ED 
scale (120 /xHz). 

For UiG, however, the symmetry of multiplets 
is broken for frequencies greater than l.SmHz ap- 
proximately, indicating that, from rotational ve- 
locities around 16kms~^, second-order effects of 
rotation become important. Nevertheless, such ef- 
fects are not strong enough to mix up multiplets of 
different mode degrees, whose centroid modes are 
separated (in average) by about 10 /iHz(0-2 multi- 
plets) and 20/LtHz (1-3) multiplets. However, the 
rotational splitting for this model is 3 /iHz and 
the TO ridges are almost two times (in average) 
more spread in the ED than those of model Uqs- 
Consequently, the identification of the mode de- 
gree is not so obvious in this case: The (1,+1) 
ridges are closer to (3,-3) as the frequency in- 
creases (around 4/iHz at iy„,i — l.SmHz, and 
2 /iHz at — 3mHz). Similarly, radial modes 
ridges are close to (2,-2) ones (less than 5 /iHz). 
This is particularly relevant for the highest fre- 
quency domain. 

But the most extreme case considered here is 
the U32 model. With such a large rotation rate, 
second-order effects including near degeneracy ef- 
fects (see Section 13. 2p become significant. From 
Figure [31 one notices that there is a blend of ridges 
which might blur the typical ED features, in par- 
ticular when few modes are detected. The scenario 
described for Uie is somehow intensified now. The 
blending is caused by ridges approaching, and even 
crossing each other in certain cases: The (1,-1) 
ridge crosses (1,0) at i'n,e — 1.5mHz. Similarly, 
(2,-2) crosses the radial modes at i^nj — l.SmHz. 
Moreover, (3,— to) are very close to any of the 
(l,m) ridges, e.g. (3,-1) with (1,+1), (3,-2) with 



(1,-1-1), and (3,-3) with both (1,+1) and (3,-^2). 
In such a case, unless one has good independent 
reasons to assume that only m=0 modes are de- 
tected, it would be extremely difficult to attribute 
a {£, m) value to some of the modes. 

6. Discussion & conclusions 

We showed that that even for mild rotation, 
second-order effects of rotation start to have signif- 
icant impact on the oscillation frequencies in the 
range of their expected maximal detection, that 
is around Vmax, i-e. the frequency at which the 
observed power spectrum reaches its maximum. 
Here, we use the empirical relati on calibrated to 
solar values (|Bedding et al.l 120041) can be written 
as ^ 

(10) 

where I'ac is the acoustic cut-off frequency defined 

by 



(11) 



where Hp is is the pressure scaleheight. We as- 
sume d the value for the solar cal ibration Vs.c,q / t'max,© 
1.70 (jBalmforth fc Gough|[l990l) . Indeed, rotation 
effects on EDs (see Figure[3]) are of the order of the 
/iHz for rotational velocities higher that 8kms^^ 
(e 0.018, /i(j/inax) ~ 0.6) for frequencies around 
J^max (see Tabled]) and higher. 

6.1. Mode identification 

From the above results, one question naturally 
arises: How can we distinguish between different 
TO ridges when they are close or even crossing each 
other? This question has no easy answer. There 
are many physical variables that must be consid- 
ered when analyzing mode excitation and detec- 
tion: geometry, visibility, intrinsic amplitudes, ex- 
citation mechanisms, etc. It is not the purpose of 
this work to go into great detail on these aspects, 
however, analysis of Fig |3] can provide some hints 
which could help with mode identification. For 
instance, notice that —to ridges spread more dras- 
tically than the +m ones, this effect being more 
significant as the rotation velocity increases. This 
would thus help to identify the m 7^ components 
of the oscillation spectra, in particular in the case 
for which this differential behaviour is stronger, 
i.e. for the £ = |to| modes. In fact, based on visi- 
bility and geometrical considerations, those modes 
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are more visible for increasing inclination angles, 
the limit case being i —>■ 90deg, that is when the 
star is being observed equator-on. Furthermore, 
it has been found empirically that observed ampli- 
tudes for p modes in main sequence stars are corre- 
lated with increasing i (Suar ez et al.,,2002i) . This 
result was then supported by theoretical simula- 
tions using non-pert urbative theory for th e oscil- 
lation computations (jLignieres et alJ 120061 ) which 
predict that the mode amplitude tends to con- 
centrate near the equator (firstly predicted by 
Clement 1998[ ). It can thus be concluded that the 



larger the rotational velocity of the star is, the 
more chances the observed ridges in EDs corre- 
spond to = and/or i = \m\ modes. 

6.2. High frequency range and surface ef- 
fects 

Second-order effects of rotation increase with 
the frequency and are therefore largest in the high 
frequency domain (above Vmax) where frequen- 
cies are sensitive to v arious surface effects (see 
ICtouDil fc DunretJl2007L for a review on these ef- 
fects). For instance, Straka et al. ( 2006[) showed 
that turbulence causes variations in the upper 
part (highest-frequencies) of EDs. In that work, 
a comparison of classical models without turbu- 
lence with models including 3D simulations of 
turbulence reveals that differences in frequency 
for radial modes of about 5 /zHz (from i/max to 
higher frequencies) between the two types of mod- 
els can exist (see EDs of Figs. 4 and 5 in Straka 
et al.). These orders of magnitude are comparable 
to those we find here for models with rotational 
velocities of about 16kms~^ in the equivalent fre- 
quency domain. 



Si milarly, it has been sh own (jTheado et al 



20051 : ICastro fc Vauclaiill2006l ) that the effects 



the model frequencies with and without includ- 
ing such diffusion processes are different for radial 
modes than for non-radial modes thereby affect- 
ing the small spacings. Such effects can reach a 
few /iHz for a model of 1.30 M©. This means 
that, rotational second-order effects shown in the 
EDs of Figure |3] can be larger than those com- 
ing from microscopic diffusion of He even for very 
small rotational velocities. Rotation can thus be 
considered as another contribution to such sur- 
face effects, like turbulence. Note however that 
rotation effects extend to smaller frequencies for 



higher rotation rates. 

Scaled small frequency spacings like ri^o (see 
Section [23} have been shown to be rather insensi- 
tive to surface effects, and therefore can probe 
properties of the i nner part, like the convec- 



tive c ore overshoot ([Roxburgh fc Vorontsovll2003 , 
20041 ). This is no longer true when rotation is 
responsible for the 'surface effects'. Indeed we 
find a similar behavior of both spacings ri^ and 
Svi^Q (Figs. [3J; ancQjD, respectively). The reason 
is that the effects of rotation on EDs are mainly 
explained by the variations in the small spacings. 
These variations become significant when the ro- 
tational velocity and oscillation frequency increase 
(see Figs. 12.1b . fc 12.1b .) but they are not compen- 
sated by the effects on large spacings which remain 
quite small. This is confirmed with a calculation of 
the ri_o removing the effect of the stellar distortion 
by the centrifugal force from the oscillations fre- 
quencies, i.e. we neglect the second-order terms in 
the perturbation treatment. Then, all the curves 
in Figure [31; collapse (Figure [IJi). 

These results imply that even at small stel- 
lar rotational velocities (from 10-15 kms~^), the 
scaled r£_„ spacings are sensitive to the stellar dis- 
tortion, therefore care must be taken when using 
them as indicators for probing the stellar interior. 
It is possible to find appropriate combinations of 
frequencies which enable to remove 'pollution' due 
to second-order effects of rotation in ri p and 5i>i 
([Pziembowski fc Goupij|[l998t [Coupil et al.[[20a 
Lochard[|200l 



6.3. An illustration: the solar-like star 

rj Bootis 

In order to illustrate how rotation affects the 
asteroseismic diagnostics of real data, we se- 
lected the well-known solar-like star -q Bootis 
(HD HR52 35) which has a measu red vsini of 
12.8 kms"^ (|Bouchv fc Carrie j[2002l ). This bright 
sub-giant GO star has been the subject of many 
observational and interpretation works in the lit- 
erature. One of the most recent works on the- 
ore tical interpretation of this star was published 
by [Carrier et al.[ (|2005[) (hereafter C05). In that 
work, a detailed asteroseismic modeling of the 
star, including rotation and atomic diffusion, was 
performed. For illustrative purposes, we com- 
puted stellar models which are approximately sim- 
ilar to some of the representative models listed in 
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Table 5 of the aforementioned paper. The physical 
parameters and the observed oscillation frequen- 
cies of the star were taken from COS. 

As shown in the previous sections, shellular ro- 
tation and uniform rotation profiles yield almost 
identical results at the scale of EDs. Thus, in or- 
der to simplify the procedure, we adopt, for the 
present illustration, the hypothesis of uniform ro- 
tation. For the convection description, the mixing- 
length formulation is used, assuming an overshoot 
of 0.2. The physical parameters and the observed 
oscillation frequencies of the star were also taken 
from C05. 

We searched for models matching the observed 
frequencies and the large spacing estimated by 
C05 ((Az/„,f) = 39.9 AiHz) with the help of EDs. 
To do so, stellar models with rotational velocities 
ranging from 12 kms^^ {i ^> 90 deg) and 48 kms^^ 
[i — > 14 deg) were computed. This range of rota- 
tional velocities ensures the use of a perturbative 
approach for the oscillation computations (see dis- 
cussion of the perturbative approach in Section H]). 

Assuming that most detected modes are m = 
modes, the best diagnostics for the large separa- 
tion {{Ai/n/) = 39.78 ^Hz) were obtained for a 
model, A4rp with similar characteristics to model 
M3 of C05, i.e., a mass of 1.70 Mq with radius R = 
2.79 R0, and an age of 2.25 Gyr (Tes = 6014 K, 
L = 9.14 L0), and surface rotational velocity of 
30kms~^ (which implies an inclination angle of 
i ~ 25°). A value a = 1.75 is taken for the MLT 
mixing length parameter. The chemical compo- 
sition is fixed at X = 0.7, Z = 0.03. Additional 
characteristics of this model are listed in Table [TJ 

The resulting ED for A^^ fits reasonably well 
the m — modes with the observations (see Fig- 
ure [5^) whereas for a lower rotation velocity set 
at 21 kms^^, the match is poorer and the stellar 
model would not be considered as a good match 
in absence of rotation (Figure \5jp), particularly 
for the £ — 1,to = modes. Even when rota- 
tion is included, the mode identification is not 
unique. Note that the (1,-1), (3,-3), and (3,-2) 
ridges are very close to (1, 0), which is also the case 
of (1, -1-1) with (3, 0). Similarly, we found that the 
ridges (2,-2) and (2,-1) are almost coincident 
with (0,0), and the (2, +2) with (2,0). Some of 
the high-frequency modes {I'n.e > 0.85 mHz) are 
better fit with m ^ ridges. For the selected 
model, these modes are hardly fit when rotation 



effects are not included in the modeling, as shown 
in panel b of Figure [S] In any case, we are aware 
that the presence of all m components makes the 
ED more dense. This increases the probability of a 
better match and non-uniqueness of the solution, 
in particular when no independent mode identifi- 
cation information is available. Model fine tuning 
should thus consider variations of both the amount 
and the shape of the rotation profile. 

Quantitatively, the effects of rotation found 
reach 5 /iHz approximately in certain cases (e.g. 
for £ — 1 modes) when considering only m = 
modes, but this is largely increased when fitting 
m 7^ ridges. In general, these effects are much 
larger than the observational uncertainties. Since 
no information on individual frequency uncertain- 
ties are given in C05, the observational uncer- 
tainties shown in Figure [S] correspond to the un- 
certainty of the observational (5i'o,2 small spac- 
ings given in that work, i.e. around 1 /iHz (^1^0,2 
corresponds to the mean separation between the 
£ — and 2 ridges). Note that it is possible that 
the present solution including rotation may not 
be unique. For instance, small variations around 
the physical characteristics of A^,, (mass, effective 
temperature, gravity, angular momentum distri- 
bution, etc.) may yield similar fits. It can be 
shown that variations in the physical characteris- 
tics of A4,, must be small in order fit the observa- 
tions as shown in Figure [5l This is coherent with 
an echelle diagram of an evolved star (in the main 
sequence) as it is predicted by model A^,,. This 
explains also that small variations of the modulo, 
around 0.01 /iHz affect significantly to the accu- 
racy of the fit. 

Depending on the star characteristics, the 
impact on precise seismic diagnostics and and 
thereby on accurate modeling varies. For very 
small rotational velocities, i.e. vsini < 5kms~"'^, 
rotation effects can be neglected. However, for 
higher vsin i, the modeling should be performed 
taking the second-order (distortion) effects into 
account. Special caution must be taken when no 
additional information about the angle of inclina- 
tion of the star is available (and/or the rotational 
splitting) , otherwise the accuracy of the modeling 
may be compromised significantly. This, com- 
bined with additional information coming from 
spectroscopy and/or multicolor photometry, may 
be helpful for the identification of modes in solar- 
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like stars (Suarez et al., work in progress using 
CoRoT data). 

This example illustrates well that for solar like 
oscillating, low and intermediate mass main se- 
quence stars with rotation rates in the upper part 
of the observed rotation rate range, the model- 
ing should be performed taking full second-order 
(distortion) effects into account. Hence, special 
caution must be taken when no additional infor- 
mation about the angle of inclination of the star 
is available (and/or the rotational splitting), oth- 
erwise the accuracy of the modeling may be sig- 
nificantly compromised if the rotation of the star 
happens to be large despite a small vsmi. 
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Fig. 3. — Echelle diagrams for the uniformly- 
rotating models Ui. Open squares correspond 
to ridges of the non- rotating stellar model A^o- 
The remaining symbols and colors have the same 
meaning as in Figure [TJ The small horizontal line 
shows the location of the i^max frequency. Stars, 
triangles, and inverted triangles represent modes 
with |m| = 1, 2, and 3, respectively. (For clarity, 
colors are used in the on-line version of the pa- 
per). 
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Fig. 1. — Perturbative (right panels) and non-perturbative (left panels) echelle diagrams computed for a 
polytropic model with four different rotational velocities (increase from top to bottom). Black and green 
symbols represent i = Q and £ = 2 modes, respectively. Blue and red symbols represent £ = 1 and ^ = 3 
modes, respectively. Filled circles represent both I — Q and m = modes. For the rest of modes, filled and 
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Fig. 2. — Small spacings (5z^o,2 (left panels) and Svos (right panels) as a function of the oscillation frequency. 
Panels on top correspond to NP oscillations. Middle panels show spacings calculated with a perturba- 
tive approach without including near-degeneracy effects. Bottom panels are similar to middle panels but 
near-degeneracy is included. Different line types correspond to spacings computed for different rotational 
velocities. Colored figures are available in the on-line version of the paper. 
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Fig. 4. — Effects of rotation on different small frequency spacings represented as function of the oscillation 
frequency. Panel (a), and (b) show the small frequency spacings 6vi^j, and Svi^q, respectively. Panel (c) and 
(d) show the scaled spacings ri o, the latter being calculated without taking second-order effects of rotation 
into account. The small vertical line indicates the location of the i^max frequency. 
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Fig. 5. — Illustration of considering rotation effects on echelle diagrams (a) for the seismic diagnostics for 
the solar-like star 77 Bootis (oscillation frequencies taken from COS) compared with the non-rotating case 
(b). The observed frequencies are represented by black filled squares. The nomenclature used to represent 
the theoretical frequencies in previous figures is also adopted here. For the sake of clarity, the to = as 
well as the radial modes are connected with lines. Note the presence of mixed modes at low frequency (near 
0.6 mHz). For the simplicity of the whole figure, these modes are not connected by lines. The uncertainty 
of the observed frequencies, 1 /iHz approximately, were roughly estimated from those of the small spacings 
given in C05 (more details in the text). This is approximately the size of the symbols. Stars, triangles, and 
inverted triangles represent modes with |m| = 1, 2, and 3, respectively. 
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